Dynamic range and sensitivity requirements of satellite 
ocean color sensors: learning from the past 


Chuanmin Hu, 1 * Lian Feng , 12 Zhongping Lee , 3 Curtiss O. Davis , 4 
Antonio Mannino , 5 Charles R. McClain , 5 and Bryan A. Franz 5 

'University of South Florida, College of Marine Science, 140 Seventh Avenue, South, 

St. Petersburg, Florida 33701, USA 

2 Wuhan University, State Key Laboratory of Information Engineering in Surveying, 

Mapping and Remote Sensing, Wuhan 430079, China 

3 University of Massachusetts, Boston, Department of Environmental, Earth and Ocean Sciences, 

Boston, Massachusetts 02125, USA 

“College of Earth, Ocean and Atmospheric Sciences, Oregon State University, Corvallis, Oregon 97331, USA 
5 NASA Goddard Space Flight Center, Greenbelt, Maryland 20771, USA 
‘Corresponding author: hu@marine.usf.edu 

Received 11 April 2012; revised 14 June 2012; accepted 10 July 2012; 
posted 23 July 2012 (Doc. ID 166494); published 24 August 2012 


Sensor design and mission planning for satellite ocean color measurements requires careful consid- 
eration of the signal dynamic range and sensitivity (specifically here signal-to-noise ratio or SNR) so 
that small changes of ocean properties (e.g., surface chlorophyll-a concentrations or Chi) can be quan- 
tified while most measurements are not saturated. Past and current sensors used different signal 
levels, formats, and conventions to specify these critical parameters, making it difficult to make 
cross-sensor comparisons or to establish standards for future sensor design. The goal of this study is 
to quantify these parameters under uniform conditions for widely used past and current sensors in 
order to provide a reference for the design of future ocean color radiometers. Using measurements 
from the Moderate Resolution Imaging Spectroradiometer onboard the Aqua satellite (MODISA) un- 
der various solar zenith angles (SZAs), typical (Atypical) and maximum (A max ) at-sensor radiances from 
the visible to the shortwave IR were determined. The Atypical values at an SZA of 45° were used as 
constraints to calculate SNRs of 10 multiband sensors at the same A,yp ica i radiance input and 2 hy- 
perspectral sensors at a similar radiance input. The calculations were based on clear-water scenes 
with an objective method of selecting pixels with minimal cross-pixel variations to assure target homo- 
geneity. Among the widely used ocean color sensors that have routine global coverage, MODISA ocean 
bands (1 km) showed 2-4 times higher SNRs than the Sea- viewing Wide Field-of-view Sensor (Sea- 
WiFS) (1 km) and comparable SNRs to the Medium Resolution Imaging Spectrometer (MERIS)-RR 
(reduced resolution, 1.2 km), leading to different levels of precision in the retrieved Chi data product. 
MERIS-FR (full resolution, 300 m) showed SNRs lower than MODISA and MERIS-RR with the gain 
in spatial resolution. SNRs of all MODISA ocean bands and SeaWiFS bands (except the SeaWiFS 
near-IR bands) exceeded those from prelaunch sensor specifications after adjusting the input radiance 
to A typica] . The tabulated A typical , A max , and SNRs of the various multiband and hyperspectral sensors 
under the same or similar radiance input provide references to compare sensor performance in pro- 
duct precision and to help design future missions such as the Geostationary Coastal and Air Pollution 
Events (GEO-CAPE) mission and the Pre-Aerosol-Clouds-Ecosystems (PACE) mission currently being 
planned by the U.S. National Aeronautics and Space Administration (NASA). © 2012 Optical Society 
of America 

OCIS codes: 280.0280, 280.4788, 010.4450. 
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1. Introduction 

In the past several decades, satellite ocean color 
radiometry has evolved from the proof-of-concept 
Coastal Zone Color Scanner (CZCS, 1978-1986) to 
modern instruments such as the Sea-viewing Wide 
Field-of-view Sensor (SeaWiFS, 1997-2010) and 
the Moderate Resolution Imaging Spectroradiometer 
(MODIS, 1999-present for Terra and 2002-present 
for Aqua). Added to this suite of optical measure- 
ments were the Medium Resolution Imaging Spec- 
trometer (MERIS, 2002-2012), the Ocean Color 
Monitor (OCM-1, 1999-present; OCM-2, 2009- 
present), and the most recent Geostationary Ocean 
Color Imager (GOCI, 2010-present). These frequent 
and synoptic observations of the global ocean pro- 
vided unprecedented opportunities to study the 
ocean’s biogeochemistry, trace ocean circulation, 
monitor harmful algal blooms (HABs) and pollutions 
such as oil spills, document coastal water quality 
changes, and help make management decisions on 
coral reef ecosystems and fisheries. 

To continue these observations without data gaps, 
future satellite ocean color missions are being 
planned by the U.S. National Aeronautics and Space 


Administration (NASA), the European Space Agency 
(ESA), and several other international agencies. 
Among the planned missions are the three tier-2 Dec- 
adal Survey missions from NASA, namely (1) Aerosol- 
Clouds-Ecosystems (ACE), (2) Geostationary Coastal 
and Air Pollution Events (GEO-CAPE) [1], and 
(3) Hyperspectral Infrared Imager (HyspIRI), each 
with ocean color capability and a unique set of science 
goals [2]. Under the auspices of the climate initiative, 
NASA is also planning for an ocean ecosystem sensor 
as part of the Pre-ACE (PACE ) mission planned for a 
2019 launch to ensure global ocean color climate data 
continuity [3] . 

Mission planning calls for careful consideration of 
the measurement and instrument requirements in or- 
der to achieve the mission goals. One of the challenges 
in the planning phase is to specify the signal-to-noise 
ratio (SNR, or sensor sensitivity) for every spectral 
band at a typical input signal (radiance) received 
by a satellite sensor. 

SNR is a critical parameter to determine the pre- 
cision of the derived ocean color data products [4] , for 
example, the surface chlorophyll-a concentration 
(Chi in mgm' 3 ). As an example, Fig. 1(a) shows 



Fig. 1. (Color online) Instrument sensitivity to resolve small changes in at-sensor radiance [(a), mW cm -2 /mV 1 sr _1 ] and surface ocean 
chlorophyll-a concentrations [(bMd), mg m -3 ] derived from several ocean color instruments, (a) CZCS at-sensor radiance at 443 nm from a 
small region (25° to 27 °N; 85 °W to 83 °W) in the eastern Gulf of Mexico (inset figure), (b)-(d) chlorophyll-a concentrations derived from 
the same region from MODISA, SeaWiFS, and MERIS-FR on 12 March 2009, respectively, using the same algorithms and soft- 
ware (SeaWiFS Data Analysis System or SeaDAS, version 6.3). The NE-SW linear line on the georeferenced MERIS image is a result 
of push broom striping noise. Note that MERIS-RR data have much higher SNRs than MERIS-FR data and therefore would not lead to the 
image noise shown in (d). 
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the at-sensor total radiance measured by CZCS in a 
2° x 2° region in the eastern Gulf of Mexico, where 
the pixelization noise (i.e., variations from pixel to 
pixel) due to quantization for this 8 bit sensor is ap- 
parent. The derived Chi images from the same box on 
12 March 2009 from three ocean color sensors re- 
vealed various noise levels [Figs. 1(b) — 1(d)] , with 
MODISA showing the lowest noise and MERIS-FR 
(full resolution) showing the highest noise [note that 
MERIS-RR (reduced resolution) data have much 
higher SNRs than MERIS-FR data and therefore 
would not lead to the image noise shown here; see 
Section 4]. Using a simple method (see Section 5), 
the RMS noises from adjacent pixels for various 
Chi values were quantified and listed in Table 1. 
For low-Chl waters (0.01 to 0.2 mgm' 3 ) the relative 
pixelization noise decreased with increasing concen- 
tration. Such noises would impact the sensitivity to 
resolve small Chi changes. For example, around 
0.1 mgm- 3 , MERIS-FR cannot resolve a 19.4% rela- 
tive change in Chi (corresponding to 0.019 mgm" 3 ), 
while MODISA may resolve changes of >4.8% 
(0.0048 mgm-' 3 ). Clearly, different SNRs from 
various ocean color sensors would lead to varying 
abilities to resolve small changes in Chi. 

In practice, it is always desirable to have high SNR 
or sensitivity to resolve small changes in the ocean. 
This is particularly important for the open ocean and 
colored dissolved organic matter (CDOM)-rich coast- 
al ocean where the water-leaving radiance is very 
small and the at-sensor signal is dominated by the 
atmosphere in all spectral bands. However, there 
is always a trade for high SNR. Due to engineering 
limitations, the dynamic range of a sensor is often 
fixed, leading to lower saturation radiance for high- 
sensitivity sensors. Figure 2 shows an example of the 
price of high sensitivity. While the red-green-blue 
(RGB) composite image of MODISA, derived from 
the low-sensitivity land bands, shows a nearly 
cloud-free image, a big portion of the ocean area is 
masked (black) in the fluorescence line height (FLH) 
[5] image derived from several ocean bands. The 
high-FLH pixels (dark green and yellow colors) were 
confirmed to contain a HAB, yet due to the saturation 
of several ocean bands, it is unknown whether the 
HAB extended into the masked area. Examination of 
the radiance data indicated that the maximum valid 
at-sensor radiance of MODISA is comparable to the 



Fig. 2. (Color online I (a) MODISA RGB composite image over 
coastal waters off Zhejiang Province of China (25°20’N to 30° 
40’N; 119°20’E to 124°E); (b) the corresponding fluorescence line 
height (FLH) image shows an HAB (dark green and yellow colors). 
Note that although the RGB image derived from the land bands 
shows minimal cloud cover, a big portion of the ocean area shows 
no FLH data due to saturation of the corresponding ocean bands. 
White lines delineate coastlines. 

“knee” radiance of SeaWiFS in the green-red bands, 
but lower than SeaWiFS “knee” radiance in the near- 
IR bands (Table 2). SeaWiFS has four detectors for 
each spectral channel; three are set to high gain 
for ocean imaging and one to low gain for land 
imaging. The “knee” occurs in the combined signal 
when the three high-gain detectors are saturated. 
Thus, thick aerosols, fog, moderate sun glint, or 
high water turbidity can cause saturation of these 
MODISA ocean bands, leading to no data in some 
important regions. This problem is more profound 
for the Japanese Global Imager (GLI [6]; also see 
Nakajima et al. [7]), as the 490 nm saturation 
radiance is only 6.4 mW cm -2 ^/m -1 sr _1 , much lower 
than the corresponding SeaWiFS knee radiance 
(8.53 mWcm-^m -1 sr _1 , Table 2). 

Thus, it is critical to design a sensor with sufficient 
SNR to resolve small changes in the surface ocean 
while not saturating the sensor under cloud-free con- 
ditions. To correct bright-target adjacency effects 
from clouds and land (i.e., the at-sensor radiance 
from the ocean is contaminated by the nearby bright 
target), it is also desirable to measure the bright- 
target (clouds) radiance for most, if not all, clouds. 
To do this we must first specify the spectral at- 
sensor radiance from typical cloud-free ocean scenes 
(^typical) anf l maximal spectral radiance from clouds 
(L max ). To cover the full range of possible conditions 


Table 1. Pixelization Noise in Chi Data Products from Three Ocean Color Instruments” 


Chi (mgm- 3 ) 

0.01 

0.02 

0.03 

0.05 

0.1 

0.2 

0.4 

MODISA 

16.5% 

11.3% 

9.3% 

7.5% 

4.8% 

2.8% 

2.8% 

SeaWiFS 

21.9% 

16.7% 

15.8% 

14.5% 

8.5% 

6.1% 

7.1% 

MERIS-FR 

74.2% 

51.5% 

40.1% 

29.5% 

19.4% 

12.1% 

8.2% 

MERIS-RR 

21.6% 

18.4% 

16.4% 

12.9% 

8.1% 

4.3% 

4.4% 


“The statistics were derived from at least >10, 000 pixels for each concentration. Chi data were derived from sensor measurements using 
near identical algorithms and processing steps defined in SeaDAS6.3. Note that MERIS data used here are both full resolution (FR, 300 m) 
and reduced resolution (RR, 1.2 km) as compared with 1.1 km of MODISA and SeaWiFS (both used Local Area Coverage, or LAC, data). 
MERIS data processed using similar algorithms but different software, such as those from the European Space Agency, showed slightly 
different RMS noise than that listed in the table. 
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Table 2. Maximum Valid Radiance of MODISA and Knee Radiance of SeaWiFS for All Ocean Bands, Determined from Measurements in 2004* 


X (nm) 

412 (412) 1 2 3 * * 6 

443 (443) 

488 (490) 

531 (510) 

547 (555) 

667 (670) 

678 

748 (765) 

869 (865) 

MODISA max valid 

19.99 

15.90 

11.56 

9.13 

6.97 

3.50 

3.31 

2.24 

1.31 

MODISA saturation* 

26.9 

19.0 

14.0 

11.1 

8.8 

4.2 

4.2 

3.5 

2.5 

SeaWiFS knee 

11.76 

10.87 

8.53 

7.29 

5.97 

3.37 

N/A 

2.46 

2.09 


“Radiances are in units of mW cm" 2 /im _1 sr -1 . MODISA maximum valid radiance is the radiance above which the linear relationship 
between the sensor-recorded digital number and calibrated radiance no longer holds. As a reference, MODIS saturation radiance refers to 
the radiance preceding sensor saturation [ 11 ] . SeaWiFS Knee radiance was derived using a second derivative analysis from the digital 
number (DN) versus radiance relationship. 

‘Wavelength (nm) in parenthesis is for SeaWiFS. 

“Data provided by Franz et al. [11]. 


from a geostationary platform from early morning to 
late afternoon, these two quantities, together with 
the SNR under typical conditions, need to be defined 
at various solar zenith angles (SZAs) in order to help 
with sensor design and mission planning (e.g., deter- 
mine integration time). 

Unfortunately, published values for SNR, L typical , 
and L max in either refereed or gray literature were 
not defined in a consistent way for previous and ex- 
isting ocean color sensors, making it difficult to learn 
from the past. Some instruments specified noise- 
equivalent radiances (NEAL) instead of SNRs [8], 
and some defined SNRs from laboratory measure- 
ments using input radiance that is not typical for 
the ocean (e.g., MODIS [9]). Even when the solar dif- 
fusers on SeaWiFS and MODIS were used to “mea- 
sure” the SNRs on orbit, the input radiances were 
different [9,10], and adjustment of the input radiance 
to realistic typical radiance over the ocean is neces- 
sary to determine the SNRs [11]. To date, for past 
and existing ocean color sensors there is a lack of con- 
sistent SNR assessment under the same radiance 
input for typical measurement conditions over the 
ocean. Such information can serve as an important 
reference to help design future ocean color missions 
such as GEO-CAPE and PACE. 

Thus, given the need for these parameters under 
uniform measurement conditions, this study was 
designed to 

(1) determine spectral L typica i and L max from 
MODISA measurements under various SZAs, 

(2) assess SNRs of commonly used ocean color 
sensors under the same L typical , and 

(3) provide tabulated references and recommen- 

dations for future ocean color missions. 

2. Data Sources 

Atypical and L max were determined from MODISA 
measurements. MODISA data were obtained from 
the NASA Goddard Space Flight Center (GSFC). 
To assure sufficient coverage under all possible con- 
ditions, MODIS data covering all of North America 
and adjacent oceans (0 °-70 °N and 135 °-45 °W) dur- 
ing different months of 2004 were selected. A total of 
255 Level-0 data granules were processed using the 
software package SeaDAS (version 6.3) to derive the 
at-sensor radiance (Level-IB), SZAs, surface Chi, 
and quality flags for each pixel (Level-2). 


SNRs were estimated from MODISA and several 
other past and current ocean color sensors, including 
SeaWiFS, MERIS-FR, MERIS-RR, CZCS, OCM, and 
GOCI. Level- 1 at-sensor radiance data covering the 
same North America ocean areas were obtained from 
various sources, including NASA GSFC (for Sea- 
WiFS LAC data, MERIS-FR, CZCS), the European 
Space Agency (MERIS-RR), the Earth Scan Lab of 
Louisiana State University (for OCM-1), and Korea 
Ocean Satellite Center (for GOCI) (Table 3). These 
data were processed in the same manner as with 
MODISA to generate at-sensor radiance, surface 
Chi, and quality flags for each pixel (Level-2). 

Although high spatial resolution sensors (e.g., 
Landsat TM and ETM+, HJ-CCD, EO-1 Hyperion) 
were designed primarily for land use, there have 
been demonstrations that they could also be used to 
study coastal blooms and water quality [12-15], 
More recently, the Hyperspectral Imager for the 
Coastal Ocean (HICO) was installed on the Japanese 
Experiment Module Exposed Facility (JEM-EF) on 
the International Space Station and has been collect- 
ing data for selected sites in the global coastal zones 
since September 2009 [16]. In addition, the geosta- 
tionary GOES Imager is equipped with a wide band 
(550—750 nm centered at 650 nm, GOES I-M Data 
Book) that could also be used for coastal ocean stud- 
ies. At-sensor radiance data from these measure- 
ments have also been obtained to estimate their 
SNRs over the ocean under the same typical condi- 
tions. Specifically, Landsat TM and ETM+ and 
EO-1 Hyperion data were obtained from the U.S. 
Geological Survey [17]. HJ-CCD data were obtained 
from the China Centre for Resources Satellite Data 
and Application (http://www.cresda.com/nl6/index 
.html). HICO data were obtained from the data ar- 
chive at Oregon State University [18]. GOES Imager 
data were obtained from the National Oceanic and 
Atmospheric Administraion (NOAA) Comprehensive 
Large Array-data Stewardship System (CLASS). The 
characteristics of these are all listed in Table 3 (see 
informations sources in [19]). 

3. /-typical and L max Determined from MODISA 
Measurements 

For each MODIS granule and spectral band, L typical 
was derived from the Level-1 at-sensor radiance from 
cloud-free ocean pixels. These pixels were defined by 


6048 APPLIED OPTICS / Vol. 51, No. 25 / 1 September 2012 


a 


m © 

w SP 


C/) 

o 


6 

w 


o'? I 2 w g 


o 

03 O 


I I 


Q 

O 

< D 

m =3 

^ M 
£ 
O 


S 4 o 

I 03 W 
00 2 50 


pq 

t“3 iH 

w ^ 


o 

I CD 
03 CD 
O | 
O lO 
03 t> 


jb 

T3 C 


,<% 
i g o 
p £ 


< 

aJ Q 

.S M 
-C H 
O « 

o 


4 


1 I < < 

g^c-o-doom^ 

5 - 5 


LO 


w 


<1 <: 


^ CO g -3 CO S ^ 

CO CD P 

4 ^ 5 


00 O 
03 t> 
i—l CD 
I I 

00 CO 
t> ^ 
03 ^ 


►» 

»0 rg 

CD 03 ^ 
^ 03 


< 

00 02 

P 


LO 
I CD 
03 00 
O | 
O 03 
03 i—l 


§ j 

>o 


£ u 
k tn 
x O 

I « 

o w 

CO 


1 

03 

03 

UJ 

CD 

CO 

LO 
03 
00 X 

4 

03 

aJ 

o 

P3 

03 

03 


r O 



CO 

i— 1 

i—l 

CD 

CO 

03 



w 


i—l LO 
O CD 
03 00 
I I 
t> 03 
03 i—l 
03 ^ 


>> 

ccS 

T3 O 
03 ’" H 
I 


1 

03 

O 

o 

°l 

LO 

o 

03 

i—l 

ccS 

03 

O 

CO 

1—1 

\ 

T3 

i—l 

03 

1—1 


o 

CO 

CO 



5 CO 


o 

o 

o 


aJ 

0 

a" 

i 

CO 


i—l 


03 

o 

1—1 

03 

| 

16 

\ 

o 

o 

ccS 

T3 

< 

o 

03 

03 

i— 1 


LO 

\ 

03 

1 





o 

i—l 





LO 






03 



& 

— 'a S 

c] i2 5 o be 

■ S n< 

CJ -S § Ss-a s. 

.2 2 f h -H 7° C 

ts tj ^ ^ .S ^ "S 
p rv O ra OJ o o 

p m £ Z; « Z; o 


1 September 2012 / Vol. 51 , No. 25 / APPLIED OPTICS 


6049 



Fig. 3. (Color online) L typ i ca i of MODISA determined from the 255 
MODIS scenes using clear water pixels under a clear sky at var- 
ious SZAs. The data are tabulated in Table 4. 


the quality control flags in the corresponding Level-2 
data. The pixels were grouped by various SZAs 
(10° ± 1°, 30° ± 1°, 45° ± 1°, and 70° ± 1°) in order 
to determine A typical as a function of SZA. For each 
group, the mean radiance from all qualified pixels 
was regarded as A typi( . al for that particular granule. 
Then, A typical for the predefined SZA was computed 
as the mean and standard deviation radiances of 
all individual A typical values from the 255 granules. 
The results are shown in Fig. 3 and listed in Table 4. 
Clearly, A typical decreases with increasing SZA. In 
reality, the at-sensor radiance is a function of not 
only SZA but also sensor angle as well as the variable 
aerosols and ocean properties. However, because our 
purpose was to determine the at-sensor radiance 
under typical measurement conditions rather than 


under exact atmospheric and oceanic conditions, 
these tabulated values (mean and standard devia- 
tions) should be able to serve as references for future 
sensors and to estimate SNRs (see Subsection 4.A ). 
Indeed, for all ocean bands in the visible, Atypical at 
SZA = 45° is significantly higher than those listed 
by NASA prelaunch sensor specifications for MODIS 
(Table 4). In contrast, the two near-infrared (NIR) 
ocean bands (748 and 869 nm) and all land bands (ex- 
cept 469 and 555 nm) showed significantly lower 
Atypical than the NASA prelaunch specifications, 
which were based on typical land reflectances. Be- 
cause Atypical at SZA = 45° was used in this study 
to determine the SNRs, such derived SNRs could 
not be compared directly with those listed in the 
NASA specifications. 

A max was derived from cloud pixels. Because all 
MODISA ocean bands (8-16) saturated over clouds, 
only land bands (1—7) were used. For each MODISA 
granule, the following steps were used to deter- 
mine A max . 

(1) Cloud pixels were selected using quality con- 
trol flags in the MODISA Level-2 data. 

(2) Of these cloud pixels, for each band, the top 5% 
and 1% brightest pixels were selected, respectively. 
The extreme outliers due to detector malfunctioning 
were discarded in this step. 

(3) The mean radiance of these qualified pixels 
was regarded as A max for the selected granule. 

Similar to the calculation of A typical , A max was com- 
puted as the mean and standard deviation of all in- 
dividual A max values from the 255 granules. These 
steps were repeated for each predefined SZA (10°, 
30°, 45° and 70°). Results for the top 5% and 1% 


Table 4. L typica | (Units: mWcm _2 //m _1 sr _1 ) of MODISA over the Ocean (0 °-70 °N and 135 °-45 °W) Determined from On-Orbit Measurements 

during All Months of 2004* 


Band # 

X (nm) 

SZA = 

10° 

SZA = 

o 

O 

CO 

SZA = 

45° 

SZA = 

= 70° 

NASA Spec 6 

Mean 

STD 

Mean 

STD 

Mean 

STD 

Mean 

STD 

8 

412 

10.27 

0.42 

9.47 

1.50 

8.07 

1.46 

4.95 

1.17 

4.49 

9 

443 

8.75 

0.42 

8.38 

1.38 

6.98 

1.33 

4.50 

1.10 

4.19 

3“ 

469“ 

7.80 

0.46 

7.43 

1.26 

6.19 

1.19 

4.01 

1.02 

3.63 

10 

488 

6.52 

0.61 

6.28 

1.10 

5.23 

1.01 

3.49 

0.87 

3.21 

11 

531 

4.36 

0.55 

4.17 

0.77 

3.55 

0.73 

2.33 

0.59 

2.79 

12 

547 

3.74 

0.69 

3.67 

0.71 

3.13 

0.66 

2.05 

0.52 

2.10 

4“ 

555“ 

3.45 

0.42 

3.35 

0.62 

2.85 

0.61 

1.85 

0.47 

2.90 

1“ 

645* 

1.72 

0.21 

1.65 

0.33 

1.39 

0.32 

0.96 

0.23 

2.18 

13 

667 

1.61 

0.21 

1.53 

0.32 

1.27 

0.30 

0.92 

0.22 

0.95 

14 

678 

1.53 

0.20 

1.43 

0.30 

1.19 

0.29 

0.87 

0.21 

0.87 

15 

748 

1.04 

0.16 

0.95 

0.22 

0.75 

0.20 

0.56 

0.13 

1.02 

2“ 

859“ 

0.55 

0.11 

0.50 

0.12 

0.40 

0.11 

0.27 

0.059 

2.47 

16 

869 

0.56 

0.11 

0.51 

0.13 

0.41 

0.12 

0.29 

0.060 

0.62 

5 C 

1240“ 

0.14 

0.042 

0.12 

0.038 

0.086 

0.038 

0.058 

0.016 

0.54 

6“ 

1640“ 

0.056 

0.019 

0.045 

0.016 

0.031 

0.016 

0.018 

0.007 

0.73 

r 

2130' 

0.015 

0.005 

0.011 

0.004 

0.008 

0.005 

0.004 

0.002 

0.10 


“The last column lists NASA specifications [9] . The mean and STD values for SZA = 45° in columns seven and eight were used as the 
typical radiance to determine sensor SNR (i.e., sensitivity) for 10 heritage instruments. 

6 The NASA prelaunch MODIS specifications were provided by Xiong et al. [9] . For SeaWiFS, the input radiances for SNR specifications 
were [9.10, 8.41, 6.56, 5.64, 4.57, 2.46, 1.61, 1.09] for the eight spectral bands at [412, 443, 490, 510, 555, 670, 765, 865 nm] (see [ 10 , 45 ]). 
“These bands were designed for land and atmospheric applications, but they could also be used to study the ocean. 
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Fig. 4. (Color online) L max of MODISA determined from the 255 MODIS scenes using cloud pixels at various SZAs. (a) L max from the top 
5% of the brightest pixels; (b) L max from the top 1% of the brightest pixels. The data are tabulated in Tables 5 and 6, respectively. 


brightest pixels are shown in Figs. 4(a) and 4(b) , 
respectively, with the L max values listed in Tables 5 
and 6 , respectively. Compared with their correspond- 
ing Ltypicai values, L max values were much higher. 
Further, there was some variability between L max 
from different images, resulting in about 10 % stan- 
dard deviation. For each band, L max decreased with 
increasing SZA. This result, together with similar 
characteristics in L typical (decreased with increasing 
SZA), may be used to help design future satellite 
missions on geostationary platforms. Note that if a 
sensor has L max values as listed in Table 5, measure- 
ments over turbid atmosphere or moderate sun glint 
such as in Fig. 2 will not be saturated. 

4. SNR at L, ypical 

A. Methods 

For a radiometer, the SNR is defined as the ratio of 
the signal electrons to the number of unwanted or 
noise electrons on the detector: 

SNR ( Slgnal/N oise ) Asi^naN^ Anise ' 

^signal = (O/hv) *t*A*rj , (1) 

where O is the power density (W nr 2 ) falling on 
the detector, hv is the photon energy (W s), t is the 
exposure time (s), A is the pixel area (m 2 ), and 77 


is the quantum efficiency of the detector at that 
wavelength. 

Three noise signals are typically identified: u photon 
or shot noise, n det or detector noise, and n rea dout or the 
noise associated with the readout electronics. These 
noise sources are not correlated, so the total noise is 

^noise S q rt Khoton + n det + Readout) ■ < 2 > 

For modern detectors and electronics, u photon domi- 
nates the noise with n 2 hoton » n 2 iet + < aflout , because 
the input signal even from a dark ocean surface still 
overwhelms the total noise [20,21]. For the purposes 
of this study, we assume that n nojse can be approxi- 
mated by n p hoton: which is proportional to the square 
root of the signal. SNR is thus a function of the 
detected signal. To properly define SNR for an ocean 
color sensor, we need to use a typical ocean spec- 
tral signal for this calculation. For this practical 
consideration, we define signal as the radiance 
(m Wcm" 2 //m 4 sr -1 ) measured from the target, and 
noise is the associated random variability [ 22 ] . This 
is different from some of the definitions in engineer- 
ing or instrument design, where u p h 0 ton collected 
from a finite field of view or falling on a detector is 
considered as the signal. In practice, all measured 
signals are digitized, so digitization round-off repre- 
sents another noise source that can cause errors in 
the derived data products [4]. In SNR evaluations, 


Table 5. L max of MODISA Land Bands over Clouds (5% Brightest Pixels) (Units: mWcm _2 pm _1 sr -1 )* 


Wavelength (nm) 

SZA = 

10° 

SZA = 

CO 

0 

0 

SZA = 

= 45° 

SZA = 

O 

O 

Mean 

STD 

Mean 

STD 

Mean 

STD 

Mean 

STD 

469 

50.15 

6.55 

46.99 

6.11 

40.45 

5.61 

23.68 

4.43 

555 

42.21 

5.50 

40.20 

5.18 

34.06 

4.86 

19.64 

4.07 

645 

37.09 

4.80 

35.09 

4.48 

30.27 

4.20 

18.02 

3.76 

859 

19.68 

1.69 

20.18 

1.28 

18.89 

1.63 

12.13 

2.36 

1240 

9.08 

0.68 

8.97 

0.74 

7.95 

0.82 

4.78 

0.97 

1640 

3.40 

0.38 

3.36 

0.37 

3.06 

0.35 

1.78 

0.38 

2130 

0.80 

0.09 

0.78 

0.13 

0.73 

0.11 

0.44 

0.09 


“Note that MODISA ocean bands all saturate over bright clouds. 
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Table 6. Same as Table 5 , but i. max was Derived from the 1% Brightest Pixels (Units: mWcm^/im -1 sr _1 )“ 


Wavelength (nm) 

SZA = 

10° 

SZA = 

CO 

O 

o 

SZA = 

45° 

SZA = 

70° 

Mean 

STD 

Mean 

STD 

Mean 

STD 

Mean 

STD 

469 

58.92 

4.28 

53.35 

6.00 

44.91 

6.13 

27.15 

4.45 

555 

50.05 

3.57 

45.59 

4.99 

38.42 

5.40 

22.49 

4.30 

645 

44.21 

3.00 

39.88 

4.36 

33.61 

4.67 

20.88 

3.94 

859 

21.50 

0.35 

21.43 

0.47 

20.48 

1.21 

14.27 

2.36 

1240 

9.90 

0.41 

9.83 

0.63 

8.67 

0.87 

5.61 

0.90 

1640 

3.78 

0.24 

3.76 

0.33 

3.35 

0.37 

1.97 

0.32 

2130 

0.93 

0.07 

0.93 

0.13 

0.81 

0.11 

0.54 

0.09 


“Note that for SZA = 10°, the listed L max values correspond to reflectances of 91%, 87%, 89%, 71%, 69%, 50, and 30%, respectively. 


this round-off-induced noise is considered together 
with inherent sensor noise without distinguishing 
the sources. Obviously, two quantities need to be de- 
termined to calculate SNR: input signal (i.e., radi- 
ance) and measured noise. 

For a given sensor, input signal affects SNR. Thus, 
to compare SNRs between different sensors, the in- 
put signal needs to be kept the same. This has been 
problematic in the past, because most published 
SNRs were referenced against different input signals 
in laboratory measurements, and Xiong et al. [9] and 
Eplee et al. [ 10 ] used onboard solar diffusers to esti- 
mate SNRs for SeaWiFS and MODIS on orbit, re- 
spectively, at predefined input signals. However, 
these predefined signals were different, making a di- 
rect SNR comparison between SeaWiFS and MODIS 
difficult. In addition, many other ocean color in- 
struments are not equipped with an onboard solar 
diffuser, suggesting that an alternative method is re- 
quired to define the input signal. Indeed, for remo- 
tely sensed imagery, the signal can be considered 
as the mean of the calibrated radiance (or gray-level 
count) from a spatially homogenous target [20,23], 
This approach was adapted in the current study to 
restrict the SNR calculations to those measurements 


where the at-sensor radiance was equal to MODISA 
Atypical a t SZA = 45° (Table 4, shaded columns). For 
sensors with different spectral bands, the tabulated 
Atypical values were interpolated. 

Similar to the restrictions on the input signal, 
noise can be represented by the standard deviation 
over a spatially homogeneous target [24,25] . The po- 
tential problem of this approach, however, is that the 
standard deviation of the selected target may contain 
real ocean or atmosphere changes within the target 
in addition to sensor noise. Thus, the target must be 
selected carefully, for example, from oligotrophic 
oceans where the horizontal gradient is small. This 
is not trivial in practice. Wettle et al. [ 26 ] demon- 
strated a method of Automated Local Convergence 
Locator (ALCL) to find a sampling size of 33 x 33 
pixels when the standard deviation approached an 
asymptotic value. Gao [ 23 ] proposed that the histo- 
gram mode of all standard deviations from prede- 
fined small windows could represent the noise. 
However, it is still unknown whether each target 
within the small window is homogenous. In this 
study, the Gao [ 23 ] approach was adapted with an 
added objective analysis of homogeneity of each 
small window. Briefly, if the maximum/minimum 




Fig. 5. (Color online) (a) SNR of the MODISA 412 nm band determined using different thresholds of maximum/minimum ratio within 
small areas (running windows on an image). The SNR values of individual images become stable (i.e., they are similar from image to 
image) when the threshold is >1.002. Below this threshold, many images did not yield any qualified small windows, and SNRs determined 
from individual images varied substantially. To minimize the potential impact of inhomogeneous water and/or atmosphere on the SNR 
determination, the max/min threshold ratio in this study was chosen to be 1.002. The same method was applied to all other bands to 
determine the threshold value to assure spatial homogeneity of the target; (b) Noise level distribution of all qualified small windows from 
a typical image. The distribution is normal (according to the Gaussian fit) when the max/min threshold equals 1.002, while it is not normal 
when the threshold equals 1.0015. 
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Fig. 6. (Color online) Distribution of noise, defined as the stan- 
dard deviation of all pixels in a given window (either 3 x 3, 5 x 5, 
or 7 x 7), as a function of window size from a simulated image. The 
image was assigned a unit value for all pixels, with added random 
noise from -1.0 to 1.0 with a Gaussian normal distribution. The 
different window sizes led to <5% changes in their histogram 
modes (numbers in parenthesis), suggesting that image noise 
(and thus SNR) is not sensitive to the choice of window size. 

radiance ratio within a small window was smaller 
than a certain threshold, the target in the window 
could be considered as homogeneous. The threshold 
was determined through an objective method 
detailed below. 

In summary, for every sensor in the study, the 
following steps were used to determine the SNR at 
input signal = A typical for each spectral band. 

(1) In each image, only clear-water pixels were 
selected to calculate statistics. For MODISA, 
SeaWiFS, CZCS, MERIS-FR, and MERIS-RR, after 
discarding all pixels associated with the quality 
control flags, the pixels with at-sensor radiance 
(L t ) within A typical ± 1 standard deviation (STD) 
and Chi < 0.07 mgm 3 were selected. A Chi value 
of 0.07 mg/m 3 is the threshold to define the most 


oligotrophic ocean waters [ 27 ] . For all other sensors 
(OCM, GOCI, TM, ETM+, HJ-CCD, GOES, HICO, 
Hyperion), clear-water pixels that meet the criteria 
of Atypical ± 1 STD were selected by visually examin- 
ing the RGB images (to exclude land and clouds) and 
the corresponding Chi images from concurrent mea- 
surements by MODISA using the same criteria as 
with MODISA. Note that the selection of clear-water 
pixels did not consider the satellite viewing angle 
even though the pixel resolution degrades from nadir 
view to side view. This is because that radiance (or 
number of collected photons) does not change with 
viewing angle for a fixed instantaneous field-of-view 
(IFOVj from a given sensor, according to the radiance 
invariance law for Lambertian surfaces. 

(2) The qualified pixels were further screened 
using the following method. All qualified pixels were 
divided into 3x3 small windows. To assure minimal 
variations within the individual windows so that the 
variations were due to sensor noise rather than due 
to real variations in the ocean or in the atmosphere, a 
threshold of maximum/minimum ratios from all 
windows was used. The windows whose maximum/ 
minimum ratios were above this threshold were dis- 
carded in the SNR calculations. STD was determined 
from each window where all 3 x 3 pixels had valid ob- 
servations, and the histogram mode of all standard 
deviations was regarded as the instrument noise 
for that particular image. SNR was calculated as 
Atypical divided by noise. Because the sensor’s SNR 
only depends on the sensor itself and should not 
change from image to image for the same A typical , by 
trial and error the threshold was determined when it 
yielded a relatively stable SNR from all images and 
when it yielded a normal-distributed histogram of all 
STDs (i.e., noise). Below this threshold, SNRs could 
not be determined from some of the individual 



Fig. 7. (Color online) SNRs of all spectral bands from the 10 multiband instruments, determined from their individual measurements at 
near-identical Atypical (right axis, dashed line). Note that the MODISA bands were divided into three groups: 1 KM (ocean bands between 
413 and 869 nm); HKM (land bands at 469, 555, 1240, 1640, and 2130 nm), and QKM (land bands at 645 and 859 nm). MERIS SNRs were 
determined from both FR and RR data. Because they were determined from multiple images, both mean and standard deviations (vertical 
bars) were plotted. These data are tabulated in Table 7. The numbers in the parentheses represent the number of images used in 
determining the SNRs. 
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images. A Gaussian fit was used to determine 
whether the noise had a normal distribution. If the 
significance level CP value) of the T-test was <0.05, 
the noise distribution was considered to be normal. 
Figure 5 shows the concept using the MODISA 
412 nm band as an example. When the threshold 
was selected as <1.002, SNRs from different images 
varied substantially and many images yielded no 
qualified small windows to start with. At 1.002, 
SNRs from all images stabilized, leading to a mean 
SNR of 1651.2 [Fig. 5(a) ] together with a normal dis- 
tribution in the STD statistics [Fig. 5(b) ] . In contrast, 
a threshold of 1.0015 would lead to a non-normal dis- 
tribution in the STD statistics. Thus, a max/min 
threshold of 1.002 was selected to screen all pixels, 
and the SNR was determined as the mean SNR from 
all examined images for this band. The procedure 
was repeated for each band and for each sensor with 
its max/min threshold ratio derived independently. 

(3) The sensitivity of SNR calculations to the vary- 
ing window size was also studied using simulations. 
A 500 x 500 artificial image with a constant pixel va- 
lue was assumed, and Gaussian-distributed random 
noise, with minimum and maximum values of -1 and 
1, respectively, was added to the image. For each win- 
dow size of 3 x 3, 5x5, and 7x7, the noise of each 
small window was determined as the STD. The his- 
togram of all STDs is shown in Fig. 6, where the his- 
togram mode appears to be insensitive to window 
size. Because under natural conditions more pixels 
selected in a small window would result in a higher 
likelihood of containing natural variability (in either 
the ocean or the atmosphere, or both), a 3 x 3 window 
should be the best choice and, thus, was chosen to 
estimate SNRs for all sensors. 

B. Results 

Such determined SNRs for all sensors are presented 
in Fig. 7 and tabulated in Table 7. Of all sensors eval- 
uated, MODISA ocean bands and MERIS-RR data 


showed the highest SNRs, with MERIS-RR SNRs 
reaching 3000 for the blue bands. The most striking 
result was that all MODISA ocean bands showed 
much higher SNRs than the NASA specifications 
(bold font in Table 7; see Xiong et al. [9]), especially 
in the blue-green bands. One reason leading to this 
discrepancy is that the SNRs in the blue-green were 
determined from Atypical much higher than the NASA 
specified input radiances (last column in Table 4). 
However, even after the input radiance from the 
NASA specifications was adjusted to Atypical (bold 
font in parentheses in Table 7), MODISA SNRs in the 
ocean bands were still significantly higher than the 
specifications. The same was true for most MODISA 
land bands; their SNRs were all significantly higher 
than the NASA specifications (after adjusting the in- 
put radiance to A ty p ica i) except for the 1640 and 
2130 nm bands. For example, at 859 nm, the NASA- 
specified SNR was 201 with an input radiance of 
2.47 mW cm -2 /im" 1 sr -1 , but our estimates showed 
an SNR of 157.4 at a much lower input radiance of 
0.40 mW cm -2 /im' 1 sr -1 . If the input radiance was 
adjusted from 2.47 to 0.40 mW cm -2 ^/m _1 sr -1 , the 
NASA SNR specification for this band would be 
80.9, much lower than 157.4 from our calculations. 
At 555 nm where A typica i was similar to the NASA 
specifications, the SNR was higher than that esti- 
mated from the solar diffuser measurements at a si- 
milar input radiance (Fig. 4 of [9]). For similar 
wavelengths, the SNRs of the land bands were much 
lower than those of the ocean bands. If these higher- 
resolution (250 and 500 m) land bands are aggre- 
gated to the same resolution as the ocean bands 
(1 km), their SNRs can increase significantly (e.g., 
compare SNRs of MERIS-FR and MERIS-RR data 
in Table 7). 

These results from MODISA estimates do not al- 
ways agree with those reported in Franz et al. [11], 
even though the A typical input radiances used in both 



Fig. 8. (Color online) Similar to Fig. 7, but these SNR values were determined from hyperspectral instruments of HICO (black symbols) 
and Hyperion (brown symbols), respectively. While the input radiance used to determine HICO SNR is similar to .^typical for MODIS (green; 
the same as Aypicai i n Fig- 7), the input radiance used to determine Hyperion SNR is higher than MODIS Aypicai i n the green, red, NIR, and 
most shortwave IR wavelengths. 
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studies are similar. While the SNRs derived here for 
the blue-green ocean bands (412-547 nm) are very 
similar to those of Franz et al. [11], the red (667 
and 678 nm) and NIR (748 and 869 nm) ocean bands 
showed 20%— 40% lower SNRs here. In contrast, all 
the land bands showed higher SNRs than those 
listed in Franz et al. [11]. We believe that this con- 
trast is because of the different reference targets 
used in the two evaluations: this study used mea- 
sured radiances from homogeneous targets as the 
reference, while Franz et al. [11] used solar-diffuser- 
based measurements as the reference. Because many 
sensors are not equipped with solar diffusers, the ap- 
proach used here should at least provide a consistent 
evaluation between different sensors. 

Next to MODISA and MERIS-RR SNRs is MERIS- 
FR in the blue-green bands, followed by SeaWiFS, 
OCM, and GOCI. For the visible bands, SeaWiFS 
SNRs all exceeded the prelaunch specifications 
(Table 7), consistent with the findings from the solar 
diffuser measurements [10], However, SeaWiFS 
SNRs in the NIR bands were lower than the sensor 
specifications or as found by Barnes et al. [28] and 
Eplee et al. [10], They were lower than OCM and 
GOCI for the same wavelengths. The disproportion- 
ally low SNRs in the SeaWiFS NIR bands have sig- 
nificant implications on atmospheric correction 
[29,30] and ocean color data products, as small noises 
in these bands would propagate to much higher 
noises in the blue-green bands through atmospheric 
correction, resulting in pixelization in Chi images [4] . 

Compared with MODISA ocean bands and 
MERIS-RR, MERIS-FR showed considerably lower 
SNRs for all bands, but most of these FR SNRs agree 
very well with those reported by Delwart [ 31 ] . This is 
understandable because MERIS-FR at 300 m resolu- 
tion collected much less photons than the MERIS-RR 
data at 1.2 km resolution. The reduced SNRs for 
MERIS FR are associated with gains in spatial re- 
solution, important for coastal zone studies. As 


expected, the proof-of-concept CZCS several decades 
ago showed significantly lower SNRs than other 
ocean color sensors. The land-oriented sensors all 
showed SNRs below 100 for all bands for the typical 
ocean radiances over nonbright targets, suggesting 
that special treatment of the data (e.g., spatial bin- 
ning) may be required when these instruments are 
used for ocean applications. 

One of the most notable results is that although 
the geostationary GOCI is at an altitude far from the 
Earth’s surface (about 36,000 km versus about 
700 km for a polar-orbiting satellite) with a higher 
spatial resolution (500 m) than SeaWiFS (1 km), 
SNRs in the visible bands were similar to those of 
SeaWiFS, and in the NIR were much higher (more 
than double) than SeaWiFS, because of the longer 
dwell times possible from a geosynchronous orbit. 
If GOCI data are binned to 1000 m resolution, the 
SNRs will be significantly higher. It is thus expected 
that GOCI will provide more precise (but not neces- 
sarily more accurate) ocean color products than 
SeaWiFS. This inference has been verified in a pre- 
liminary study using only several images [ 32 ] but 
will be tested thoroughly once GOCI data become 
routinely available. 

SNRs of the hyperspectral instruments of HICO 
and Hyperion were also estimated in the same man- 
ner. Because of the limited data availability, the SNR 
calculations were restricted to input radiance as 
close to MODISA L typical as possible, but not exactly 
the same as with the above calculations. Figure 8 
shows the SNRs of HICO and Hyperion at the corre- 
sponding input radiances. Note that the HICO input 
radiance was very similar to MODISA L typica ] for all 
bands, while the Hyperion was designed for land 
imaging and the input radiance was much higher 
than MODISA Atypical for wavelengths >600 nm. For 
wavelengths <600 nm, HICO SNRs were >200 and 
reached >400 at 500 nm for the given input radiance. 
This result agrees well with the sensor specifications 



Fig. 9. (Color online) (a) SNR at different radiance ( L t ) inputs (solid curves), modeled from SNR at one L t input using a square root 
function. Also plotted are the SNR values derived at different (A,) from direct MODISA measurements (discrete symbols); (b) MODISA 
and SeaWiFS SNRs determined from measurements over the ocean (solid symbols) at input radiance of MODISA Atypical at SZA = 45°. 
These SNRs are compared with those from the sensor specifications after adjustment of the input radiance to the same Atypical (empty 
symbols). Also shown in the figure are the MODISA/SeaWiFS SNR ratios (stars) for the same A typ i ca i input radiance. 
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[ 16 ] and is significantly higher than the Hyperion 
SNR even if the input radiance of Hyperion is much 
higher. The result suggests that while HICO can 
be used for ocean studies at original spectral and spa- 
tial resolutions, spectral and/or spatial binning of 
Hyperion may be required to improve SNRs for the 
dark ocean. 

Except for the hyperspectral sensors, SNRs of 
all sensors were evaluated in this study at near- 
identical input radiance (MODISA A typical at 
SZA = 45°) so their performance could be compared 
with each other directly. On the other hand, to 
compare them with published SNRs of other sensors 
at different input radiances, some adjustments 
are necessary, because SNR depends on the input 
signal. 

F igure 9(a) shows an example using three MODISA 
ocean bands at 443, 547, and 678 nm. As expected, all 
SNRs increased with increasing input radiance. Re- 
sults for other bands were similar. These sensors 
are designed to be shot-noise limited and therefore 
the SNR should be proportional to the square root 
of the input signal. Using this simple model, we can 
predict the SNR at another input radiance using 
the SNR at A typical (MODISA, SZA = 45°). Figure 9(a) 
shows how the predicted SNRs (solid curves) agree 
with those determined from measurements (symbols) 
for the three MODISA ocean bands. Comparison of 
the predicted and measured SNRs showed <10% dif- 
ference for most MODISA bands, with better perfor- 
mance in the visible ocean bands than in the NIR 
ocean bands and the land bands, possibly due to 
the higher SNRs in the former. The result suggests 
that SNR at an arbitrary input radiance, when it is 
within ±50% of Atypical* can be predicted from the 
known SNR at A typical , or vice versa, with uncertain- 
ties generally <10%. 

This result is also useful when comparing SNRs 
obtained in this study with other published SNRs 
at different input radiances. For example, after the 
SNRs of MODISA and SeaWiFS from the sensor 
specifications were adjusted to force the input radi- 
ance to change from the specifications (Table 4) to 
MODISA A typical at SZA = 45° using a square root 
function of radiance [Fig. 9(b), empty symbols], they 
could be compared directly with those determined 
from MODISA and SeaWiFS measurements at 
Atypical [Fig. 9(b), solid symbols]. The NASA SNR spe- 
cifications after adjustment of the input radiance are 
listed in Table 7 (bold font in parenthesis). For 
MODISA, all ocean bands except 678 nm showed 
significantly higher SNRs than the sensor specifica- 
tions (35% — 160%). Likewise, all SeaWiFS visible 
bands (except 555 nm) showed significantly higher 
SNRs than the sensor specifications (12% — 67%), 
but the NIR bands showed much lower (~30%) SNRs 
than the specifications. Overall, under the same 
Atypical input radiance, MODISA SNRs are 2. 1-4.7 
times the SeaWiFS SNRs for the visible bands, 
and about 4.4 times the SeaWiFS SNRs for the 
NIR bands. 


It has been shown that the performance of the Sea- 
WiFS NIR bands has degraded over the mission span 
(1997-2010) [10]. However, evaluation of SeaWiFS 
SNRs using measurements from 1998, 2004, and 
2010 resulted in virtually the same SNRs, with their 
long-term differences smaller than the standard 
deviations. Thus, the SNRs presented in Table 7 
may represent time -independent values for future 
references. 

5. Impact of SNR on Ocean Color Data Products 

The sensor sensitivity to detect small changes in 
ocean properties is a direct function of the sensor’s 
SNR. Consistent with the decreasing SNR trend 
from MODISA to SeaWiFS, the corresponding Chi 
data products showed increasing pixelization noise 
(Fig. 1). While the noise could be visualized from 
the images in Fig. 1, it was quantified using a similar 
scheme to that for SNR estimates [ 33 ] . For a number 
of images, noise was estimated from all valid pixels 
(after discarding those associated with various qual- 
ity control flags) whose Chi values equaled the 
predefined Chi (0.01, 0.02, ..., 0.4 mgm -3 ). For each 
3x3 window surrounding the qualified pixel, RMS 
error was estimated. Then, for all qualifying 



Fig. 10. (Color online) MERIS-FR chlorophyll-a concentrations 
on 12 March 2009 from a small region in the eastern Gulf of Mexico 
derived from (a) the default band-ratio OC3 algorithm and (b) a 
new color index algorithm [33]. Inputs to both algorithms were 
the spectral R IS data derived using the atmospheric correction ap- 
proach in SeaDAS6.3. Note that this is approximately the same 
region as shown in Fig. 1_, but without map projection, so the ver- 
tical striping on the right edge is apparent. Note the significant 
reduction in noise from (a) to (b), where eddy circulation features 
are clearly revealed. 
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windows, the mean of the RMS errors was regarded 
as a result of the pixelization noise. The ratio be- 
tween the mean RMS error and the predefined Chi 
value was used to represent the relative pixelization 
noise (in percentage). Results are listed in Table 1. 
For each predefined Chi, relative noise increased 
from MODISA to SeaWiFS to MERIS-FR. Relative 
noises for MERIS-RR Chi were comparable to those 
of MODISA. For each sensor, relative noise decreased 
with increasing Chi and became stable around 
0.2— 0.4 mgm' 3 . 

Such SNR-related Chi noise has been explained by 
Hu et al. [4] using both simulations and SeaWiFS ob- 
servations. Because the atmospheric correction used 
the NIR bands to determine aerosol type and thick- 
ness and then extrapolated to the visible bands, a 
small error in the NIR, induced by sensor noise or 
digitization round-off noise, could propagate to the 
visible bands and become much larger in the result- 
ing water-leaving radiance (L w ) or remote sensing 
reflectance (R IS ). The R rs errors were further propa- 
gated in the various products through the bio-optical 
inversion algorithms. For the blue/green band ratio 
Chi algorithm [ 34 ] that was used in SeaDAS as 
the default, relative errors in the blue/green ratio de- 
pend on the magnitude of the ratio. For oligotrophic 
oceans with Chi <0.1 mgm 3 , the blue/green ratio 
could reach 8.0. Thus, a small error in the denomina- 
tor (green band) would lead to a much larger relative 
error in the ratio, resulting in larger relative errors 
in the derived Chi. Other data products (e.g., absorp- 
tion and backscattering coefficients, attenuation 
coefficient) may respond differently to the sensor 
noise, yet the principle is the same, i.e., low sensor 
SNR will lead to higher noise in the data products. 

As demonstrated in the Chi data product for the 
noise-induced errors, most errors originated from 
the NIR atmospheric correction bands rather than 
from the blue-green bands that were used in the 
bio-optical inversion [4] . This explains why MODISA 
Chi showed much lower pixelization noise than Sea- 
WiFS Chi (Table 1), as MODISA SNRs in the NIR 
bands (748 and 869 nm) are much higher than Sea- 
WiFS SNRs in the corresponding NIR bands (765 
and 865 nm). However, the same principle cannot ex- 
plain why MERIS-FR Chi showed much higher pix- 
elization noise than SeaWiFS Chi or why MERIS-RR 
Chi showed much higher pixelization noise than 
MODISA Chi. Because MERIS NIR bands at 754 
and 865 nm have similar SNRs to SeaWiFS SNRs 
(for MERIS-FR data) and to MODISA SNRs (for 
MERIS-RR data), one would expect similar perfor- 
mance in the Chi data product derived using the 
same algorithms, but the results here did not support 
this expectation. MERIS Chi data obtained from the 
ESA using similar atmospheric correction and bio- 
optical inversion algorithms but different software 
(BEAM) showed similar results. It was speculated 
that this abnormal behavior might be related to the 
push broom design or the offset stability of MERIS. 
MERIS has an active offset compensation loop (OCL) 


to modulate the “effective zero” signal differently 
across the spectral bands, which will accentuate 
along-track modulations during atmospheric correc- 
tion. However, we compared SNRs and Chi pixeliza- 
tion noise estimated from MERIS-RR data during 
2010 (OCL-ON) and 2005 (OCL-OFF) and found no 
difference. Currently the relatively high pixelization 
noise of MERIS Chi data from such a high-SNR sen- 
sor remains a puzzle to be further investigated. 

Because the Chi pixelization noise mainly origi- 
nated from R IS errors due to imperfect atmospheric 
correction, once an improved atmospheric correction 
algorithm is developed to reduce the R IS errors or a 
bio-optical inversion algorithm insensitive to the 
R rs errors is developed, noise in the product may be 
significantly reduced. Recently, a band-subtraction 
color index (Cl) algorithm was proposed for Chi retrie- 
vals over oligotrophic oceans (Chi < 0.25 mgm -3 ), 
and the algorithm was proven to be much more toler- 
ant to R ts errors induced by the imperfect atmo- 
spheric correction [33]. Figure 10 shows an example 
of the improved MERIS-FR Chi product using such 
a band-subtraction algorithm. Note that the R rs in- 
puts to the band-ratio algorithm [Fig. 10(a) ] and the 
new band-subtraction algorithm [Fig. 10(b) ] were 
identical, yet the noise in the default band-ratio 



C*hl ( ms nr-') 



Fig. 11. (Color online) Sensitivity of MODISA normalized FLH 
(nFLH) in detecting Chi changes at low concentrations. Data were 
extracted from a 10° x 10° (60 °-50 °W, 30 °-40 °N) region near the 
North Atlantic Gyre. nFLH was derived from normalized water- 
leaving radiance at 667, 678, and 748 nm. (a) nFLH appears 
to monotonously increase with increasing Chi for Chi > 
0.08 mgm -3 (winter) and for Chi > 0.14 mgm -3 (summer). nFLH 
values are much lower during summer due to nonphotochemical 
quenching; (b) The sensitivity of nFLH to resolve small Chi 
changes is dependent on the number of pixels used in calculating 
the averages. The annotated numbers indicate the number of pix- 
els used to calculate the mean nFLH values. It appears that the 
relationship between nFLH and Chi for December 2011 became 
stabilized for N > 500. 
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Chi has been reduced significantly using the new al- 
gorithm. The figure here used MERIS-FR data to 
show the impact of algorithms on product noise for 
the same input data. For inter-sensor comparison, 
the same atmospheric approach implemented in Sea- 
DAS was used to process all ocean color data. If a 
different approach such as the artificial neural net- 
work (ANN) designed for MERIS [ 35 ] or other semi- 
analytical inversion algorithms [36—39] were used, 
the resulting noise in the Chi products could be signif- 
icantly lower. 

For sensors with a finite dynamic signal range, 
there is always a trade between SNR (sensitivity) 
and saturation radiance. It is typical that the SNR 
for one or two bands drives the sensor design. In the 
MODIS case the high SNR required for the FLH 
bands drove the MODIS design; the telescope and 
other optical components were scaled to achieve the 
desired SNRs for those bands and, in doing so, all 
other bands had excess SNRs (e.g., significantly 
higher than 1000 for the blue-green bands; see 
Table 7). To maximize the SNR for open ocean scenes, 
the saturation level was lowered. This often resulted 
in signal saturation in the FLH and atmospheric 
correction bands over river plumes, intense algal 
blooms, or thick aerosols (e.g., Fig. 2), providing no 
usable data for those regions. Even so, at the current 
sensor sensitivity MODISA FLH could not resolve 
Chi fluorescence changes at low Chi concentrations 
from individual images. 

An example is shown in Fig. 11 for the North 
Atlantic. For December 2011, even after averaging 
>10,000 Level-2 1 km pixels, normalized FLH 
(nFLH) could only resolve Chi down to 0.08 mgm' 3 . 
Averaging more pixels, especially after extensive bin- 
ning in both space and time such as in Behrenfeld 
et al. [40], may further resolve Chi to about 
0.05-0.06 mgm -3 , yet the large nFLH variance for 
a given Chi value [STD bars in Fig. 11(a) ] indicated 
that on individual images the ability of nFLH to re- 
solve low Chi variations was significantly degraded. 
This was confirmed by comparison between indivi- 
dual nFLH and Chi images (not shown here) and 
by the sensitivity test in Fig. 11(b) , where the nFLH 
versus Chi relationship became stabilized only after 
averaging >500 pixels. This ability further decreases 
during summer due to nonphotochemical quenching 
of the fluorescence efficiency (black symbols in 
Fig. 11(a) ; also see Behrenfeld et al. [40]). The current 
MODISA SNRs appeared sufficient to derive nFLH 
for studying phytoplankton physiology in the global 
ocean at large spatial-temporal scales, because many 
data points could be averaged to increase its sensitiv- 
ity [40], yet for individual images at original resolu- 
tion (1 km), increasing the SNRs in the FLH bands 
while maintaining a sufficient signal dynamic range 
for a variety of water types may be a technical 
challenge for future ocean color sensors. Such cap- 
ability is particularly important for detecting and 
quantifying blooms in optically complex waters 
[41]. Indeed, the noise level in the MODIS nFLH 


data, as measured by the standard deviation term in 
Fig. 11(a) for the December 2011 case, is about 
0.003— 0.004 mW cm -2 /on -1 sr -1 , corresponding to a 
Chi change of ~0.1 mgm -3 [i.e., from 0.1 to 
0.2 mgm 3 in Fig. 11(a) ]. Thus, on individual images, 
current MODIS SNRs cannot resolve nFLH changes 
corresponding to Chi changes of 0. 1 mg m -3 when the 
default FLH algorithm [5] is used. An alternative al- 
gorithm based on radiative transfer theory [ 42 ] may 
improve such nFLH sensitivity in detecting Chi 
changes. 

In summary, low sensor SNRs would lead to higher 
data product noise. The atmospheric correction 
bands in the NIR play a dominant role in propagat- 
ing sensor noise to the visible bands and data pro- 
ducts through atmospheric correction [4,43], and 
thus should be emphasized in sensor design. Finally, 
noise in the current data products can be reduced 
with improved algorithms even when the sensor’s 
SNRs remain the same. 

6. Discussion 

A. Estimation of SNR from Satellite Measurements 

One objective of the SNR assessment in this study 
was to put all sensors under the same reference input 
radiance from realistic measurements so they could 
be compared directly. The approach used in calculat- 
ing the SNRs, as detailed in Subsection 4.A, was 
adapted from a community-accepted method but ad- 
justed using an objective method to restrict to ocean 
targets as homogeneous as possible through the use 
of max/min ratio thresholds. This is certainly not the 
only way to calculate SNR. An alternative approach 
could use extrapolation rather than max/min ratio. 
To evaluate whether different approaches could re- 
sult in similar results as shown in Table 7, an extra- 
polation-based method was implemented and tested. 
In this method, the criteria to select clear-water pix- 
els (Chi < 0.07 mgm 3 ) under L typical were identical 
to those used in Subsection 4. A . For each of the 20 
MODISA images, SNRs for a band were calculated 
for different window sizes ranging between 2x2 
and 7x7. The mean SNR of the 20 images for a win- 
dow size was regarded as the SNR for that window 
size. An extrapolation of the mean, as a function of 
window size, to the window size lxl was regarded 
as the sensor SNR. Such extrapolation-derived 
MODISA SNRs are comparable to those shown in 
Table 7 to within 20%, suggesting that both methods 
could be used to determine sensor SNRs. However, 
the standard deviations of the mean SNRs (i.e., 
variability of SNRs across individual images) from 
this alternative approach were significantly higher 
than the method (Fig. 5) used to derive Table 7, in- 
dicating that the alternative method depends heavily 
on the selection of the image. Thus, the general con- 
sistency between the two independent methods and 
the significantly lower cross-image variability in the 
former method suggest that the results listed in 
Table 7 could be used as valid references to compare 
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sensor specifications and performances, and the 
same method could be used for other or future sensor 
evaluations to ensure consistency. 

B. Other Factors Affecting SNRs 

Instrument engineering teams have found various 
approaches to meeting SNR specifications, some of 
which can introduce data artifacts that are not in- 
cluded in the SNR evaluations here. For instance, 
MODIS used 10 detectors along the track to ensure 
sufficient dwell time to achieve the desired SNRs, 
i.e., slow scan rates. This introduced cross-track 
striping due to residual calibration errors between 
detectors. Similarly, push broom sensors such as 
MERIS employed multiple detectors across the 
track, which resulted in along-track striping errors. 
In the SNR calculations, these errors were excluded 
by using the max/min threshold approach (Fig. 5), for 
otherwise realistic variations in the ocean and/or at- 
mosphere would have been regarded as noise. The 
striping noise represents a significant error source 
for many sensors (MODIS, MERIS, Landsat, etc.) 
that has not been addressed in this paper because 
it is primarily a cross-detector calibration problem. 
SeaWiFS avoided this problem by using four detec- 
tors in a time-delay-integration scheme but had to 
spin the telescope at a faster rate. Its bilinear gain, 
separated at the knee value, was implemented to 
account for both high sensitivity at low radiance 
and high saturation levels for land and atmosphere 
applications. 

The assessment did not attempt to differentiate 
whether the noise originated from the instrument- 
inherent noise [Eqs. (1) and (2)] or from digitization 
round-off noise, with the latter determined by the 
number of digitization bits and the maximum 
radiance. In practice, to maximize sensor perfor- 
mance, one digital number is often set to be equiva- 
lent to the instrument-inherent noise or NEAL, 
because there is no reason to digitize noise or keep 
the noise way below the digitization round-off error. 
Thus, the results presented here represent the im- 
age-based evaluations of the SNRs of many sensors 
under similar radiance inputs over relatively homo- 
geneous waters without differentiating the noise 
sources. 

C. Implications for GEO-CAPE, PACE, and Other Future 
Missions 

The ultimate goal of this study was to use past and 
current ocean color sensors as references to help 
design future ocean color missions including GEO- 
CAPE and PACE. From the results and discussions 
above, substantial experience and lessons have been 
learned from the past, including the following. 

(1) MODISA L typical values at SZA = 45° are often 
different from previous sensor specifications. They 
also change with different SZAs. 

(2) MODISA is much more sensitive (higher 

SNRs) than SeaWiFS, yet the low saturation radi- 


ance in the red and NIR bands sometimes resulted 
in saturation even over cloud-free ocean scenes. 

(3) Both MODISA and SeaWiFS showed higher 
SNRs than prelaunch specifications for most spectral 
bands. 

(4) It is critical to have sufficient SNRs in the 
atmospheric correction bands. In contrast, the SNR 
requirements on the visible bands can be relaxed. 

(5) Noise in the current data products is primarily 
due to the algorithmic approach. Improved algo- 
rithms can lead to significantly reduced product 
noise and significantly enhanced product precision. 

GEO-CAPE is a geostationary mission with recom- 
mended hyperspectral sensor specifications span- 
ning from the UV (350 nm) to the shortwave IR 
(2130 nm). One of the advantages of geostationary 
sensors is their ability to start to achieve the integra- 
tion time required to meet the SNR requirements. A 
primary measurement goal of GEO-CAPE is to quan- 
tify short-term changes in ocean properties from 
early morning to late afternoon when SZA is <70°. 
The L typical values at various SZAs, as listed in 
Table 4, can serve as references for GEO-CAPE 
SNR specifications. Likewise, the L max values listed 
in Tables 5 and 6 at various SZAs can serve as refer- 
ences to determine the saturation radiance and sen- 
sor integration time under various measurement 
scenarios. Although Table 6 is preferred, the differ- 
ence between L max for the 5% and 1% brightest pixels 
provides a reference for the trade space between sen- 
sitivity and saturation. Note that these L typical and 
L ma x values were determined using MODISA multi- 
band measurements. For hyperspectral sensors, 
Ltypical and L max at other wavelengths can be derived 
from the tables using interpolation and extrapolation 
and using the extraterrestrial solar irradiance as a 
reference. 

It is clear that MODISA SNRs are sufficient to re- 
solve very small changes (<10%) even at extremely 
low Chi concentrations (~0.03 mgirT 1 2 ; Table 1). 
When the new band-subtraction algorithm was used, 
the sensitivity was even higher to resolve changes of 
<5% for Chi between 0.01 and 0.4 mgm~ 3 ([33]; 
Fig. R)). Thus, MODISA SNRs can be used as a tem- 
plate to define GEO-CAPE specifications. Because it 
is important to be able to detect small changes dur- 
ing early morning and late afternoon at SZA = 70°, 
SNRs should be defined at this SZA. 

Based on the MODISA results shown here, and 
based on the radiative transfer simulations of Wang 
[43], we believe that the following SNRs should pro- 
vide sufficient sensitivity for the GEO-CAPE mis- 
sion: >1000 for X between 350 and 720 nm, >600 
for 720-900 nm, and >100—200 for the shortwave in- 
frared (SWIR) bands. For all spectral bands, satura- 
tion levels should be set to or interpolated from those 
listed in Tables 5 and 6, which will avoid saturation 
over land and all but the brightest clouds. The 
simulations of Wang [ 43 ] considered two cases for 
maritime and tropospheric aerosols with an aerosol 
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optical thickness at 865 nm of up to 0.2, where the 
SNR requirement was estimated to be slightly lower 
than listed above. Note that these recommended 
SNRs are much lower in the visible and slightly low- 
er in the NIR, but significantly higher in the SWIR 
when compared with MODISA SNRs determined 
from ^typical a t SZA = 45° (Table 7). Because most 
noise in the data products results from the atmo- 
spheric correction bands in the NIR and SWIR as 
opposed to the blue-green bands [4,43], the perfor- 
mance of GEO-CAPE in resolving small changes in 
ocean properties should be at least comparable to 
that of MODISA if the same atmospheric and bio- 
optical inversion approaches are used with the 
NIR bands. If SWIR bands are used for atmospheric 
correction of both MODISA and GEO-CAPE, signifi- 
cantly improved performance should be achieved 
from GEO-CAPE over MODISA because of the in- 
creased SNRs in the SWIR bands. 

Of particular interest for GEO-CAPE is the ability 
to derive changes in chlorophyll-a fluorescence quan- 
tum yield with repeated diurnal measurements. Can 
GEO-CAPE with similar SNRs to MODISA resolve 
such changes? Morrison [ 44 ] showed that when the 
surface photosynthetically available radiation (PAR; 
/.i mole photons m~ 2 s -1 ) changed from 1000 to 600 
(corresponding to SZA = 60° and 70°, respectively), 
the fluorescence quantum yield nearly doubled (Fig. 9 
of [44]). For waters with Chi = 0.15 mgm -3 , this 
quantum yield change would result in nFLH changes 
of ~0.004 mWcm^/zm^sr' 1 [Fig. 11(a) ], distin- 
guishable from the background noise level [0.003- 
0.004; see standard deviation in Fig. 11(a) ]. This 
estimate is based on individual images. If the Chi 
is higher, or if many pixels are averaged for these 
low-concentration waters, smaller changes in fluor- 
escence quantum yield can also be quantified 
with the current MODISA SNRs and the default 
algorithms. 

7. Conclusion 

The trade space between sensor SNR (or sensitivity) 
and dynamic range (or saturation radiance) needs to 
be determined carefully in order to quantify small 
changes of ocean properties (e.g., Chi) under most 
measurement scenarios. Using measurements of 
MODISA and other ocean color sensors together with 
a new objective method, typical (Atypical) and maxi- 
mum (A max ) at-sensor radiances were determined 
under different SZAs, with SNRs determined from 
uniform ocean targets at the same input radiance 
as A ty pi ca i. These results make it possible for the first 
time to have direct cross-sensor comparisons of SNRs 
under the same measurement conditions. The tabu- 
lated results also serve as references to compare with 
other existing sensors and to help design future sen- 
sors such as the geostationary GEO-CAPE or the 
polar-orbiting PACE. 
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